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Abstract 

The basic mechanism underlying the generation of Moire' patterns in 
line scanner data acquisition systems is examined. A general expression 
is developed in terms of typical system parameters for the reproduced 
image of such systems and the interaction of the image spectrum; the 
raster frequency and digital sampling frequency of the A/D conversion 
process are discussed and examples given. System design requirements 
for avoiding Moire' pattern generation and two-dimensional aliasing are 


discussed. 
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MOIRe' PATTERNS AND TWO-DIMENSIONAL ALIASING IN 
LINE SCANNER DATA ACQUISITION SYSTEMS 


I. Introduction: 

The flying spot scanner and optical line scanner are widely used 
for converting visual or other two-dimensional radiant signals to 
electrical form for further processing. When this processing is 
accomplished by means of a digital computer, A/D conversion of the 
scanner analog output signal is required. Thus in effect, the original 
two-dimensional spatial signal is sampled along two approximately 
orthogonal axes: by the scanner raster in a direction normal to the 

scan lines, and by the A/D sampling in a direction parallel to the scan 
lines. Figure 1 shows a typical system of this type used for remote 
sensing of the earth's surface from an aircraft. The transverse motion 
of the field of view of the sensor is produced by a rotating mirror, 
while the forward motion is produced by the translation of the platform 
carrying the sensor. A/D conversion may be accomplished at the time of 
measurement or later at a data processing center from an analog 
recording of the scanner signal. 

The measured signal can be reconstructed as a sampled raster by 
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Figure 1 



Basic Line Scanner Configuration Used in Remote Sensing of 
Earth's Surface. 
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intensity modulating a scanning light beam and recording the signal on 
photographic film or viewing It on a digital display. Figure 2 shows a 
typical picture recorded and reproduced In this manner. The signal was 
recorded by a line scanner operating In the Infrared portion of the 
spectrum (0.62 to 0.66 |jm). The scanner used to make the picture was 
carried In an aircraft flying at an altitude of 215 meters over typical 
farmland In central Indiana. It is evident in this picture that there 
has been distortion introduced in the form of Moire 7 patterns through 
much of the left center of the picture. Since the analog signal is low 
pass filtered to half the sampling rate of the A/D conversion process, 
the resulting Moire* patterns arise principally from the fact that 
sampling and reconstruction were carried out on a signal having signif- 
icant spatial frequency components In excess of half the scan line or 
raster spatial sampling frequency. Thus, these Moire 7 patterns, because 
of the low pass filtering prior to the A/D conversion process, are 
primarily a form of one-dimensional aliasing. Figure 3 shows the magni- 
tude of the two-dimensional Fourier transform of Fig. 2. Aliasing along 
the y-spatial frequency axis, corresponding to aliasing about integer 
multiples of the raster sampling frequency, is clearly evident. If the 
low pass filtering is inadequate, then the resulting Moire 7 patterns 
arise from two-dimensional aliasing. 

This type of distortion is particularly insidious because it cannot 
be removed a posteriori without serious loss of information content of 
the picture. Its presence indicates that serious undersampling of the 
image has occurred and that it is impossible to accurately reproduce the 
picture by interpolating between sample points to obtain a continuous 
image. Limitations of this kind are particularly important in connection 



Figure 2 
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Figure 3. Modulus of the Two-Dimensional Fourier Transform of Figure 2. 
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with the problems of correcting geometrical distortion, registration of 
data from images recorded at different times, and combining geometrical 
and spectral characteristics of images to improve classification 
accuracy. The mathematical model of the process whereby this type of 
distortion is produced and certain of the methods for handling it are 
discussed in the following sections. 

II. One -Dimensional Aliasing 

The one dimensional aliasing problem will be considered first as an 
introduction to a more precise two-dimensional aliasing problem formula- 
tion. The generation of Moire* patterns by the analyzing of patterns 
containing periodic intensity variations has been known for centuries 
[12] and much analysis has been performed on this subject [4, 8, 10, 11]. 
This process can be modelled mathematically as the multiplication of 
two or more two-dimensional intensity functions containing spatial 
frequency components that are nearly equal to one another. In the case 
of a line scanner (or a TV system) where no A/D processing is involved, 
the resulting image is taken as the product of the intensity function of 
the image and the scanner intensity function. For this case a Moire* 
pattern will result whenever there are discrete spatial frequency compo- 
nents in the image that approach the spatial frequency components of 
the scanner function. 

The basic mechanism whereby Moire* patterns are -generated is most 
easily seen by considering the two-dimensional frequency spectra of the 
image and scanner functions. Let 0(x,y) be the intensity variation of 
the image as a function of the spatial coordinates x and y, cp(x,y) be 
the intensity variation of the scanner function, and y(x,y) be the 
resulting intensity function of the scanned image. These functions can 



- 7 - 


be related by 


y(x,y) - 0(x,y) cp(x,y) 


( 1 ) 


Taking the two- dimensional Fourier transform of Eq. 1 gives the fre- 
quency spectrum of the scanned Image as a two-dimensional convolution 
of the transforms of <?(x,y) and cp(x,y); i.e., 


r(f x ,f y ) » **(f x ,f y ) **0U X . f y ) , ( 2 ) 

where f^ and f are the spatial frequency variables of the Fourier 
transform. In order to more clearly illustrate the mechanism by which 
the new frequencies (Moire 7 patterns ) are generated, it is convenient to 
consider the image function to be infinite in extent and the scanned 
image function to be the product of a raster function, cp / (x,y), that is 
infinite in extent and a window function, w(x,y), that is unity over the 
viewing area of the scanner and zero elsewhere. With these stipulations 
Eq. 1 can be rewritten as, 

y(x,y) » 6 (x,y) cp'(x,y) w(x,y) 

- 3 _1 {^ (f x »f y ) ** 0'(f x , f y )} w(x,y). (3) 


Thus the scanned image is just the inverse transform of the convolution 
of the image with the infinite extent raster function multiplied by the 
window function to limit the extent of the final image. 

In order to illustrate in a simple manner the mathematical basis of 
the mechanism whereby Moire' patterns are generated, assume for the 
moment that the image intensity function and the line scanner function 
both vary in amplitude in a sinusoidal manner. The frequencies and 
spatial orientations for these two functions need not be the same. The 
two-dimensional Fourier transform of this type of function consists of 
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a pair of impulses corresponding to each sinusoidal component. Two 

representative functions and their transforms are shown in Fig. 4 . 

The product of the raster function and the image function has a 

transform given by the double convolution of the individual transforms. 

This is readily seen by inspection (since only impulses are involved) to 

consist of four impulses as shown in Fig. 5. The two impulses near the 

origin represent a low-frequency sinusoidal component. This component 

arises from one-dimensional aliasing of the y-spatial frequency 

components of the image intensity function about the y-frequency 

components of the raster function, and correspond to what is normally 

considered a Moire 7 pattern. This Moire 7 pattern has a frequency of 

1 

f - ( f 2 + f „ 2 - 2 f . f cos a ) 2 (4 ) 

m ' i r i r ' ' 

and an angle with respect to the y-axis of 

1 f, sin a 

P = sin — , (5) 

m 

and would appear as shown in Fig. 5. If the high frequency component, 
f h in Fig. 5, is not filtered out in the reconstruction process, it will 
also be visible as a sinusoidal intensity variation in a direction 
parallel to a line connecting the two impulses at f^ shown in Fig. 5. 
Further study of Fig. 5 makes it clear that there is considerable 
Interdependence between the Moire 7 pattern and both the frequency and 
orientation of the generating signals. For example a change in orien- 
tation of one of the generating patterns can lead to a change in both 
the frequency and orientation of the Moire 7 pattern. 

With this brief introduction a more precise formulation of the two- 
dimensional aliasing problem will now be made. 
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III. Scanner- Reproducer System 

The block diagram of a typical line scanner and associated 
reproducing system is shown in Fig. 6. The variables in the system 
of Fig. 6 are defined as follows: 

x^,y^ - coordinates of scanned image 

8 (*l # y ) - intensity function of scanned image 

g (xi.yi) - aperture function of scanner 

u^t) - analog time varying output signal of scanner 

h(t) - impulse response of signal conditioning and/or low pass 
filter of A/D converter 

u 2 (t) - analog output signal of signal conditioner 

T^ - period of sampling interval in A/D converter 

x,y - coordinates of reproduced image 

g r (x,y) - aperture function of reproducer 

© r (x,y) - intensity function of reproduced image 

The output signal from the scanner at any instant of time is the 
integral of the product of the aperture function and the image intensity 
function and is given by 

N 

2 

00 00 

u l(t) » £ J J e 8 ( x 1 *y 1 ) . g 8 [*i“ v (t- nT ). y 1 “ nAyJdXj dy x (6) 

<*qd • oo 

N 

no- 

where 

H + 1 ■ number of scan lines in scanned image 
v ® velocity of scanner aperture in x^ direction 
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Ay » distance between centers of consecutive scan lines 
measured in y^ direction 

T = time duration of one scan line. 

The output of the signal conditioner, Ug(t), is the convolution of 

u^ (t ) with h (t ) and can be expressed as 


) ■ Uj_(t) * h(t) 


N 

2 


J 1 IJ S « %,[<•*!- At y^niy] 


N 

“”’2 


h (t ) dx^ dy^ dT 


(7) 


The output of the sampler of the A/D conversion process, u^(t), can 


be approximated by 


A 

2 


u (t) « u 2 (t) £ dCt-offj) 


( 8 ) 


A 

o-- 

2 


where 


A+ 1 « number of samples of Ug(t) t a k®n during one scan line. 

The reproduced image is given by a running average of the repro- 
ducer aperture function and the reproducer input signal, u^(t), as 

M 

2 

e r (x»y) - £ J u 3 (0 8 r [ x- v ( c “ mT )» y - ®^y] dt (9) 

M 

m«i- r ----- 

2 

where m+1 *> number of scan lines in reproduced image. Substituting 
Eq. 7 and 8 into 9» interchanging the order of Integration and 
summation and noting that for synchronous reproduction of the raster 


/l All integrals are to be interpreted as having infinite limits 
unless otherwise specified. 



the indices m and n are identical gives 


N A 

2 2 


e r < x ’ y) ° 1 I JJJ 9 . ( V y l ) S s [v y(0( V T '” T) ’ y l' n “'] 


N A 

n “-2 “-2 


g^x- vfcrt^ - nT), y-nAy] h(t) dXj dy 1 dt . 

( 10 ) 

Considerable insight into the significance of the various factors 
in Eq. 10 can be gained by transforming to the spatial frequency domain. 
Taking the two-dimensional Fourier transform of Eq. 10 gives 

N A 

t £ -j27l[f v(orr -nT) + f nAyJ 

f ,f ) - G (f ,f ) 7 V e x 1 y 

r' x’ y ' r' x* y' L 

N A 

m*. (x i' y i ) *.[ x i ~ v(a v t - ni >- y i - niy ] 


h(r) dXj dy^ dr . 


(ID 


The triple Integral of Eq. 11 can be recognized as a triple convolution. 
Rewriting Eq. 11 using convolutional notation 

N A 

£ 2 -j 27T[f v(cfT - nT) + f nAy] 

A'V * G r< f A> l l • 1 

N A 
n--- a~- 

[{ 0 s [v(aT 1 - nT), nAy] ** g g [-v(aT 1 - nT),-nAy]} 

* h(aT 1 - nT)J (12) 

where "**" denotes a two-dimensional convolution. Equation 12 can be 
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further simplified by writing the double summation as a double Integral 
over a sum of impulses as follows 


(f ,f ) 
r x' y' 


G r (f x' f y ) 1 J [{ e s [vx ‘ y] ** ^- vx ’- y ^} * h(x) ] * 



A 

2 

^ 5 (x - + nT , y-nAy) 



-j 27f(vf x + fy) 
e y dx dy 

(13) 


and using the fundamental properties of convolution and multiplication 
of the Fourier transform to give 


- G r<W [»<V> * Q < f x' f y>] 

(14a) 


where 


N A 

2 2 


Q(f x ,f y ) - JJ £ £ 6(x-0flT 1 +nT, y-nAy) e ' j 2TT ( vf x * f f y y )dxdy 


N A 
2 °"'2 


sin7l(N + l) T*f s in 7T (A + 1 ) T. vf sin7l(N + l) Ay f 


sin 7TT vf 


sin 77 T. vf 
1 x 


sin 71 Ay f 


(14b) 


In order to gain further insight into Eq. (l4a), consider the 
limiting case in which N-*oo and A-*eo. For this case Eq. (l4b) becomes 


u. q ( f x .fy)- i I 5(f x -q- v )^ l ^y-p - '<«>' 


A-* 00 


Q^-oo 


n«>oo 


Using Eq. 15, 


Eq. 14a can be written as 
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^4 (f ,f ) 
r' x’ y 


G (f ,f ) 
r N x* v 7 

vTj Ay 


I I H [ v(f x*V ) ]' 9, . (f x-l“-v. 

C£»-od n*”00 


f -f- ) 

y Ay 


x _ T^v ^ * ’^ f y“Ay^] * 

From Eq. 16 It is evident from the double sunnations that the 
reproduced image contains a superposition o£ frequency translates, at 
intervals of ^ and ~ along the f^ and f^ spatial frequency axes, 
respectively, of the product of the transfer functions of the analog 
signal conditioner and/or low pass filter of the A/D processor, the 
scanned image, and the scanner aperture. Each of these translates is 

> 

weighted by the reproduction aperture transfer function. Figure 7 shows 

how the magnitude of *0* r (f x ,f^) in Eq. 16 might look. The degree to 

which portions of the replicated spectra overlap into the spectral 

band - < f < ^7— , - < f < — is the degree to which 

2Ay ~ y — 2Ay 2T.jV - x -2T v 

aliasing occurs and the degree to which image distortion occurs. If 
there are discrete components in the spectrum H(vf x ) o (f x , ) 

G (-f ,-f ) near any harmonics, of the raster frequency or, ~ 

of the A/D sampling process then these will appear in the reproduced 
image as low frequency components and correspond to what are usually 
called Moire' patterns. The absence of Moire' patterns in a reproduced 
image does not necessarily imply that no distortion is present, but only 
that it is not present as a spurious narrow-band isolated component. 

Since it is impossible to remove aliasing from a sampled image, it 
is important to consider how to avoid the occurrence of this undesirable 
phenomenon. Such avoidance requires that the spectrum 




Figure 7. Partial Two-Dimensional Replication of the Original Image Spectral 
Components which Contribute to the Spectrum of the Reproduced Image. 
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H(vf x ) ^ g (f x ,f y ) G fl (-f x ,-f y ) have no f x components above^;, and f y 

components above 7 ^-. The factor *0* 8 (f x ,f y ), the spectrum of the 

original linage. Is generally not under the control of the experimenter. 

Therefore, control of the spectrum must be accomplished by means of the 

temporal transfer function of the signal conditioner and/or low pass 

filter of the A/D process, H(vf x ), and the scanner aperture transfer 

function, G (f ,f ). By causing the product of H(vf ) G (-f ,-f ) to 
s x y x 8 x y 

act as a low pass filter having a cutoff frequency of^~ in the f 

1 ^ v 

direction and — r— in the f direction, all aliasing will be eliminated. 

Achieving an approximation to the ideal low pass filter character- 
istic in the f x direction is somewhat more easily accomplished because 
the spectral transmission in this direction is a function of the product 

H(vf ) G (-f ,-f ), while in the f direction all the filtering must be 
x s x y y 

provided by G (-f ,-f ). The primary constraints governing the choice 
8 x y 

of the impulse response of the analog signal conditioner (the A/D low 

pass filter), h(t), are that it be cauBal and decay rapidly with time 

having no significant secondary lobes so as not to introduce any ghost 

images in the x-direction. However, an additional constraint is imposed 

upon the scanner aperture (or point spread function) for physical 

realizability, namely that it must not be negative. To simplify the 

discussion which follows, it will be assumed that G (f ,f ) is symmetric 

8 x y 

with respect to the f^ and f axes. Thus the problem of choosing 

G s (f x ,f y ) so that the necessary low pass filtering in the f direction 

is obtained will be considered first. 

Although not optimal, an excellent choice for g (x,y) from the 

standpoint of realizability and shape of its frequency spectrum, 

G (f ,f ), is a Hamming window function, 

8 x y 
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2 

g 8 (x*y) ■ 0. 54 + 0.46 cos ~(x 2 +y 2 ) , for x 2 +y 2 £ -jp (17) 

o 

■ 0 , otherwise 


where z q is the diameter of the circular aperture. This function, 

g (o,y), and its Fourier transform, G (o,f ) are shown in Fig. 8. The 

8 8 y 

2 

first null of G (o,f ) is at f « and this can be used to specify the 
8 y y z o 

width of the aperture z . Thus placing the first null of G (o,f ) at 

o s y 

specifies that z q ■ 4Ay, or that the scanner aperture be four times 

the width of the spacing between the centers of adjacent scans of 

the image. Having specified G (f ,f ), H(vf ) can be chosen so that 

s x y x 

the product G (-f ,-f ) H(vf ) approximates a low pass filter in the f 
s x y x x 

direction with a cutoff frequency of - t—-. Very often it is desirable to 

2T X V 

have equal spatial resolution along both axes, for such cases T^v = Ay. 

To illustrate the effectiveness of this filtering action in 
eliminating Moire' patterns, an image was generated, assuming t^v = Ay, 
having two distinct sinusoidal frequency components. Figure 9 s was 
generated by a low frequency sinusoid 

f L (x,y) «* cos 0.03 Ay (x sin 1° + y cosl°) 

and Fig. Sfo is the Moire 7 pattern produced by a high frequency sinusoid 


f H (x,y) = cos 0.986 Ay x . 

Figure 9c is the sum of f L (x,y) and f H (x,y) after sampling and recon- 
struction using an Impulse approximation for the scanner aperture 
g (x,y). Figure 9 d is the sum of f. (x,y) and f„(x,y) after sampling 

8 L H 

and reconstruction using the Hamming function of Eq. 17 for g (x,y). 

s 

It is clear that the low frequency Moire 7 distortion sinusoid (i.e., the 
vertical band of Fig. 9c) is almost completely removed in Fig. 94* 





Figure 9a. Low Frequency Sinusoid. Figure 9b. Moire Pattern Produced 

3 by a High Frequency Sinusoid and 

the Scanner Raster. 



Figure 9c. Sum of the Low Frequency Figure 9d. Output of a Hamming Spatial 

Sinusoid of Figure 9a and the Moird Filter with the Components of Figure 9c 

Pattern of Figure 9b. as Input. 


Figure 9. Example of Moire Pattern Removal by a Hamming Spatial Filter. 
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IV. Conclusions 

Moire 7 patterns appearing in the output of line scanner and 
digitally sampled data acquisition systems are a form of two-dimensional 
aliasing resulting from sampling at frequencies less than twice the 
highest spatial frequency components present in the original data. The 
existence of such patterns indicates that aliasing has occurred but 
their absence does not indicate an absence of aliasing, only that there 
are no discrete spatial frequency components in the data near a harmonic 
of the raster or sampling frequency. It is not possible to remove the 
distortion due to aliasing although specific components such as 
individual Moire 7 patterns can be removed by narrow band spatial fre- 
quency filters. The only adequate way to handle the distortion due to 
two-dimensional aliasing is to avoid it by proper design of the scanner 
aperture and by low pass filtering of the data before A/D conversion 


is carried out. 
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